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Simulating the vacuum properties of QCD with dynamical quarks 
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The vacuum properties of lattice QCD with staggered quarks are investigated by an efficient simulation method. 
I present data for the quark condensate with flavor number Nf = 0, 1, 2, 3, 4 and many quark masses, including 
the vacuum energy in the chiral limit. Obvious sea quark effects are observed in some parameter space. I also 
describe a mechanism to understand this and a formula relating the chiral condensate and zero modes. 



1. INTRODUCTION 

Lattice QCD is a fundamental theory of non- 
perturbative strong interactions, which should 
tell us everything about low energy hadronic 
physics. Here I would study the vacuum prop- 
erties of full QCD (unquenched theory). As is 
known, the most important issues are confine- 
ment, chiral-symmetry breaking and sea quark 
effects [|. For the pure gluonic sector, there 
have been active studies of the topology of the 
vacuum H. For the fermionic sector, the inter- 
play between zero modes fjl and chiral-symmetry 
breaking should be further investigated. 

Earlier studies of this issue employed quenched 
approximation, where sea quarks configurations 
are completely ignored. Quenched Chiral Per- 
turbation Theory (QCPTh) predicts existence of 
some possible extra chiral In terms , which 
do not exist in full QCD. If this is the case, the 
quenched theory is wrong and will lead to misbe- 
havior of some quantities. 

However, it is difficult in practice to visualize 
this misbehavior, due to the ambiguity in chiral 
extrapolation. To calculate < ipip > \m=o, a mass 
term rmpTp is conventionally introduced, and then 
extrapolation at some small set of m has to be im- 
plemented using some unjustified fitting function. 

In §, the Hybrid Monte Carlo (HMC) algo- 
rithm was used to simulate unquenched config- 
urations with 4 flavor Kogut-Susskind quarks. 
This is a big step towards the full theory. In 
most simulations, however, sea quarks (for the 



configurations) are usually distinguished from the 
valence quarks (for the propagators). By fixing 
m sea , and varying m va i an ce, the formula 

(i>ip(m va i ence , m sea . V)) 



J[dU}e- s °detA{m sea )J2 



3V 2m, 



A 2 +m 2 



J[dU]e-^detA(m sea ) 

is used and then extrapolated to m va i ence = 0. 
Unfortunately, this is not theoretically consistent. 

In the following sections, I would describe a 
possibly better solution to these problems. 

2. ALGORITHM 

HMC is not efficient for investigating the chi- 
ral properties. There exists an efficient algorithm 
for investigating the vacuum properties of lattice 
gauge theory with dynamical fermions, named 
the microcanonical fermionic average method |Q. 
This algorithm is based on the definition of an 
effective fermionic action 



S? ff (E,m,N f ) 



In 



J[dU][detA(m, U)] N !/ 4 8{S g {U) - 6VE) 



J[dU]S(S g (U)-6VE) 

where E is the gauge energy. The thermodynam- 
ical properties can be obtained from the deriva- 
tives of the partition function 



(1) 



Z= \ dEe- s 'ff( E < m ' N f^. 



(2) 
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Here S e ff is the full effective action related to 
S^ff, E and density of states n(E) by 



S Bff (E,m,N f ,0) = -\nn(E) 



6/3VE + S? ff (E,m,N f ). 



(3) 



This method has been tested in the Schwinger 
model H, applied to QED in 3 dimensions ||, 
QCD in 3 dimensions 0, QED in 4 dimensions 
Q , and generalized to the fermion-gauge-scalar 
models ||. For m ^ 0, it has been cross-checked 
with HMC and the data are in good agreement. 
(Fairly speaking, it is more convenient to use 
HMC to calculate the spectrum and matrix el- 
ements in the non-chiral regime) . One prominent 
advantage is that this algorithm might be still 
valid in the chiral limit, ft can also be carried 
out on small computers or workstations. 

The results to be presented will be the first 
application to QCD in 4 dimensions. 

3. RESULTS 

I describe only the 6 4 lattice data, while those 
on larger lattices will be given elsewhere. 

n(E) in (^) can be directly evaluated by nu- 
merically integrating the quenched SU(3) data: 



\nn{E) 
V 



= 6 



/ dE'/3{E',N f = 0) + const., (4) 
Jo 



as shown in Fig . [j] Fig. |2| only shows the ef- 
fective fermionic action as a function of E. The 
shapes of the curves for other Nf are quite sim- 
ilar, but the scales are quite different. We have 
approximately S^jj oc Nf . Using the saddle point 
technique, the mean plaquette energy < E > 
as a function of (5 is found, as shown in Fig. 
H for m = and different Nf. To my knowl- 
edge, < E > m= o m QCD has never been ob- 
tained before. As seen in this figure, < E > in- 
creases with Nf. The sea quark effects become 
most important around (3 rj 5.5, which corre- 
sponds to < E >« 1.5. We can understand this 
phenomenon by looking at Fig. ^, from which 
one can observe maximum slope around E « 1.5. 
The lighter the quark mass, the bigger the slope 
is. Therefore, we have a mechanism: the sea 



quark effects reach their maximum where the ef- 
fective fermionic action versus the gauge energy 
has a maximum slope. In Fig. |], the results for 
< > at m = 0.05 for different Nf are pre- 
sented. Again, one sees < •tptfj > decreases with 
Nf and the sea quark effects become dominant 
in the same region. Detailed analysis will be pre- 
sented soon. 

4. OUTLOOKS 

The estimate for < ipip > | m= o remains to be 
done. An improved method [[[o| (also 11 1) uses 
l^ijj) — 37rlimA^o For a finite V, the eigen- 
value density p(0) = due the absence of exact 
zero mode. Again, some assumed fitting function 
has to be used. 

In ]T3), Azcoiti, Laliena and I discovered that 
the chiral condensate can be directly obtained 
by analyzing its probability distribution function, 
constructed from the eigenvalues. We derived 



7,, ELi 1 

< iplp > |m=0 = y < y > 



(5) 



Here ai is the ith zero of the Bessel function Jo. 
Only a few smallest eigenvalues are necessary for 
this calculation. It is worth mentioning that no 
A or to extrapolation should be done. The only 
thing to do is the finite size analysis. One can also 
understand apparently the relation between zero 
modes and chiral-symmetry breaking. If chiral 
symmetry is spontaneously broken, the eigenval- 
ues Xi should behave as ai/V for large volume. 
We have checked [U it in the Schwinger model, 
where exact results are available. I do believe the 
method is promising for QCD. 
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Figure 2. —S^ff as a function of E for Nj = 2. 



Figure 4. < ■0V' > as a function of /3 for m = 0.05. 
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